Formation of collective spins in frustrated clusters 
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Using magnetization, specific heat and neutron scattering measurements, as well as exact calcula- 
tions on realistic models, the magnetic properties of the La3Cu2V09 compound are characterized on 
a wide temperature range. At high temperature, this oxide is well described by strongly correlated 
atomic 5=1/2 spins while decreasing the temperature it switches to a set of weakly interacting and 
randomly distributed entangled pseudo spins 5* = 1/2 and 5 = 0. These pseudo-spins are built over 
frustrated clusters, similar to the kagome building block, at the vertices of a triangular superlattice, 
the geometrical frustration intervening then at diflferent scales. 



I. INTRODUCTION 

A lot of work has been devoted recently to the under- 
standing of the peculiar magnetic behaviour of extended 
networks of spins on triangle or tetrahedron based lat- 
tices. These elemental configurations of spins, especially 
when they are corner-sharing, imply strong geometrical 
frustration for antiferromagnetically coupled Heisenberg 
spins, i.e. inability to simultaneously minimize the pair 
interactions. The case of 5*= 1/2 spins is especially ap- 
pealing because it should enhance quantum physical be- 
haviour. The archetypical frustrated network of spins in 
2 dimensions is the kagome lattice, which is expected to 
stabilize, at low temperature, a spin liquid state One 
of the most striking feature of this state is the presence 
of many low lying singlet states 3- This feature could 
be interpreted in a short-range resonating valence bond 
picture of the kagome ground state 0] . Few experimental 
realizations of kagome lattice with 3=1/2 spins are avail- 
able. Among these, the recently studied hebertsmithite 
Q does not present any magnetic transition down to the 
lowest temperature but its fluctuating ground state is 
still under debate since it does not meet the theoreti- 
cal predictions 0, 0, S @] ■ The theoretical studies have 
been undertaken on ideal systems but exact solutions can 
only be obtained for a limited number of spins (up to 36) 
[1, 7] . This is one of the reasons why it should be inter- 
esting to study experimentally a system constituted of 
isolated frustrated clusters of spins, for which the exact 
calculations are fully relevant. 

The physics of spin clusters has attracted much atten- 
tion by itself. Quantum dynamics of mesoscopic magnets 
and decoherence effects by specific environments for in- 
stance can be investigated in organo-metallic molecular 
magnets with metal clusters stabilizing a collective spin 
magnetically isolated by the organic ligands. A number 
of investigations thus followed the discovery of a stair- 
case magnetization hysteresis loop associated with the 



quantum tunneling of collective spins trough a magnetic 
anisotropy energy barrier d, Q. Clusters with small 
total spin and large Neel vector also raised strong in- 
terest. This is the case of the polyoxovanadate, short- 
named Vi5, constituted of 15 ions. Its magnetic 
properties, showing quantum processes under sweeping 
field [13], are accounted for by the formation of a collec- 
tive pseudo-spin 5=1/2 from the three weakly coupled 
5=1/2 spins on the central triangle of the cluster [ll|. 
The study of spin clusters within oxide compounds was 
undertaken more recently. At variance with molecular 
compounds, the magnetic screening then is much less ef- 
fective, which leads to non-negligible inter-cluster inter- 
actions and therefore allows the study of the coupling of 
the collective entities. For instance, in the oxide com- 
pound La4Cu3MoOi2, a trimer of 5=1/2 spins is stabi- 
lized in each triangle of Cu^^ at the vertices of a square 
lattice [H, [1^ Il4l The inter-cluster couplings in 

La4Cu3MoOi2 is two orders of magnitude smaller than 
the intra-cluster ones, and lead to the onset of an antifer- 
romagnetic long-range order at 2.6 K. In both V15 and 
La4Cu3MoOi2 compounds, the splitting of the 2 funda- 
mental 5=1/2 spins doublets might be attributed to a 
departure from the trigonal symmetry. A last remark- 
able oxide compound whose properties are dominated by 
spin clusters is Na2V307 _16, 17, 18, 19). The clusters 
are built on rings of 9 antiferromagnetically interacting 
V^^ ions of 5=1/2 spins, piled up in tubes with inter- 
ring ferromagnetic interactions. One 5=1/2 spin out of 9 
is left active in the low temperature paramagnetic regime 
through a dimerization process of the other spins of the 
ring, characterized by a broad range of V-V interactions. 
The proximity of a quantum critical point at H=Q is sug- 
gested from the complex magnetic behaviour observed 
with respect to the magnetic field H at low temperature 
T. 

The layered oxide La3Cu2V09 [13, [2l[, parent of the 
La4Cu3MoOi2 compound, provides a model at the in- 



tersection of geometrical frustration and spin cluster re- 
search areas. In the magnetic layers, the Cu^+ ions form 
planar clusters of antiferromagnetically coupled 8 and 
9 5*= 1/2 spins, arranged on 4 corner-sharing triangles. 
The geometrically frustrated 9-spin cluster is very simi- 
lar to the building block of a kagome network. As will be 
shown in the following, a collective state of spin is con- 
structed from the S=l/2 spins wave functions of each 
cluster while decreasing the temperature. The result- 
ing collective pseudo-spins involve a greater number of 
individual spins than in the V15 or La4Cu3MoOi2 com- 
pounds. Like in La4Cu3MoOi2, the inter-cluster cou- 
plings in La3Cu2V09 are two orders of magnitude smaller 
than the intra-cluster ones, and their influence is trace- 
able in the experimental low temperature range. More- 
over, the clusters themselves are arranged on a triangular 
lattice, which gives the opportunity to study the influence 
of geometrical frustration at different spatial and energy 
scales in the same compound. 

In this article, following a preliminary report [12], we 
will focus on the individual behavior of these pseudo- 
spins, their interplay when inter-cluster interactions be- 
come effective will be the subject of a forthcoming paper 
[2^. The synthesis and crystallography of La3Cu2V09 
polycrystalline sample will first be described, with special 
emphasis on the Cu/V substitution. Then, the magnetic 
measurements will be presented, in a wide temperature 
range, following the construction of the collective pseudo- 
spin, and also in high magnetic fields. The next two sec- 
tions will respectively present the results of specific heat 
measurements and of inelastic neutron scattering. The 
last section before concluding will be devoted to a dis- 
cussion of these experimental results in the light of exact 
calculations performed on realistic models of spins clus- 
ter. 



II. EXPERIMENTAL 
A. La3Cu2V09 sample characterization 

1. Synthesis and structure 

A polycrystalline La3Cu2V09 sample was synthesized 
by a sol-gel method. The stoichiometric metallic cations 
were dissolved in nitric acid before being complexed by 
addition of EDTA in a controlled pH solution. This so- 
lution was polymerized and then heated at 700 ° C to 
eliminate the organic constituents. The resulting powder 
was annealed during 15 days at 1010 ° C. The sample was 
first checked to be a single phase using X-ray diffraction. 

Neutron powder diffractograms were then recorded on 
the high resolution powder diffractometer D2B at the In- 
stitut Laue Langevin with a wavelength of 1.59 A using 
a Ge monochromator (cf. Fig. [2]). The absence of struc- 
tural phase transition was checked from 3 to 300 K. The 
crystallographic description of the sample, obtained from 
the Rietvelt analysis of the neutron diffractograms with 
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FIG. 1: (Color online) (a) La3Cu2V09 structure projected 
along the b axis: LaOs/a layers alternating with (Cu/V)03/3 
ones with the O polyhedron drawn around Cu and V 
atoms, (b) La3Cu2V09 structure projected along the c axis: 
(Cu/V)03/3 layers showing the 9 Cu^"*" clusters at the vertices 
of a triangular lattice. La, V, Cu and O atoms are represented 
by large green, yellow, blue and small red circles, respectively. 



the FuUprof program [2j| , was found consistent with the 
results of Ref. the La3Cu2V09 oxide compound 

crystallizes in the hexagonal P63/m space group with 
lattice parameters a=b=14. 395(2) A, c=10. 656(2) A at 
300 K. As shown in Fig. [T] (a), La06/3 layers alternate 
with (Cu/V)03/3 layers underlying the 2D structure of 
the compound. In these, the Cu^+ ions, the only mag- 
netic ones with a S'=l/2 spin, are localized on three in- 
equivalent sites Cu(2) Cu(3) and Cu(4) (following the 
site labeling of ref. \2u ) having different coordination en- 
vironments. These are forming 9-spin planar clusters, 
constituted with 4 corner-sharing triangles. The clusters 
are centered at the vertices of a 2D triangular lattice (cf. 
Fig. [1] (b)). Within the planes, 0-coordinated V atoms 
arc intercalated between the clusters. 

A Cu/V substitution, of vanadium ions for cop- 
per Cu'^+ ions, dictated by the global charge neutrality 
requirement that fixes the stoichiometric proportions of 
the different ions, occurs in the very small ratio of 
Notice nevertheless that this in average corresponds to 
one Cu/V substitution over ^ of the 9-spin planar clus- 
ters. We will detail the nature and consequence of this 
substitution below. 
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FIG. 2: La3Cu2V09 powder neutron diffractogramm 
recorded at 3 K on D2B, the corresponding refined pattern 
and their difference. The refined structural parameters are 
presented in table H] 



Cu(3) 




In non substituted clusters, the neighbor Cu-Cu dis- 
tances within the cluster range between 3.35 A and 3.8 
A, whereas the shortest in-plane distances between the 
peripheral Cu of the neighboring clusters is larger and 
equal to 4.55 A. In comparison, the distance between 
two adjacent (Cu/V)03/3 layers is 5.3 A. 

Within the (Cu/V)03/3 layers, the first neighbor 
5=1/2 spins of each cluster are interacting through su- 
perexchange via an intermediate O. Concerning the inter- 
cluster magnetic interaction within one layer or between 
adjacent layers, possible longer superexchange paths can 
be effective through one or two O. 



2. Cu/ V Substitution 

In order to understand the magnetic properties of the 
compound, it was found essential to take into account 
the Cu/V substitution and the related modifications of 
the cluster structure. 

To get insight into the distribution of the V in the clus- 
ters, refinement of the powder neutron diffraction data 
at 3 K was made allowing V substitution on each cop- 
per sites, alternatively. Although the fit is slightly better 
when the Cu(2) site (central triangle) is substituted in- 
stead of Cu(3) and Cu(4) sites, in agreement with previ- 
ous results 21] (see Table U and Fig. further analysis 
was estimated necessary to unambiguously characterize 
the Cu/V substitution. 

Using the VASP package [l^, Density functional the- 
ory (DFT) methods applied to electronic structure calcu- 
lation, based on the local-density approximation (LDA) 
and Projector Augmented Wave (PAW) method, were 
used for this purpose. The DFT calculations were per- 
formed with periodic boundary conditions on one unit 
cell containing 130 atoms and including 2 clusters located 
in the two (Cu/V)03/3 layers at z=l/4 and z=3/4. 

The non substituted (and hence non stoichiometric) 



FIG. 3: (Color online) DFT calculated positions of the Cu 
(blue circles), V (large yellow circles), and O (small red cir- 
cles) atoms in the 9-spin clusters (a) and 8-spin clusters (b). 
The labeling of the different superexchange interactions is 
used in section UlI Al 



structure was first optimized, and found very close to the 
one determined from powder neutron diffraction. Differ- 
ent calculations were then performed replacing, in the 
optimized structure, one Cu by a V, alternatively on 
each Cu(2), Cu(3) and Cu(4) site of the cluster. The 
large energy difference between these calculated struc- 
tures first confirmed the localization of the substitution 
on the Cu(2) sites : Ecn(3) - EcuH) > -Ecu(2) by ~ 2 eV, 
with £'cu(i) the energy of the structure with the Cu/V 
substitution on the Cu(i) site. 

In addition, the calculated structure reveals large de- 
formations around the V substitute. The two first neigh- 
bor O of the V get closer to it and further away from the 
neighboring Cu (cf. Fig. [3|). A consequence of this is the 
apparent greater isotropic displacement parameter B^so 
of the neighboring oxygen on site 0(1) in the Rietveld 
analysis of the neutron diffractogramms (see Table U. 

Calculations with two V/Cu substitutions on the same 
cluster were finally carried out and compared to those 
with only one substitution. The obtained energy differ- 
ence between these two structures, of the order of 1 eV, 
is in favor of only 1 Cu/V substitution per cluster. The 
resulting reasonable assumption of a negligible number 
of 7 (or less) Cu clusters therefore leads to a microscopic 
population of 2/3 of 9-Cu clusters and 1/3 of 8-Cu/l-V 
clusters. 
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TABLE I; Refined structural parameters from Rietveld analysis of neutron powder diffraction data of La3Cu2V09 at 3 K. Biso 
is the isotropic displacement parameter in A^. The agreement factors are 3.38 for the Bragg R-factor and 2.27 for the Rf-Bragg 
factor. 
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FIG. 4: La3Cu2V09 magnetic susceptibility from M/H mea- 
surements: in the high temperature range at 6 T using the 
BS magnetometer and at 0.1 T using the Quantum Design 
MPMS magnetometer; in the low temperature range using 
the purpose built SQUID magnetometer in 0.01, 0.1, 0.2 T 
(open circles), in good agreement with the higher tempera- 
tures measurements (plain circles). A discontinuity around 2 
K is pointed out by the vertical line. 



B. Magnetic properties 

1. Magnetic susceptibihty 

Magnetization measurements were performed at the 
Institut Neel, on two purpose built magnetometers (BS) 
using the axial extraction method in the temperature 
range from 1.6 K to 800 K and magnetic field up to 
10 T and on a more sensitive commercial Quantum De- 
sign SQUID magnetometer (MPMS) in the temperature 
range 2-350 K and magnetic field up to 5 T. Additional 
measurements at lower temperatures, from 0.23 K to 4 K, 
using a dilution inserted device were obtained on another 
purpose built highly sensitive SQUID magnetometer. 

The thermal variation of the initial magnetic suscepti- 
bility X is displayed in Fig. |4l x deduced from the 
M vs. H measurements in small magnetic fields such that 
the isothermal magnetization M varies linearly with the 
magnetic field H. An anomaly is observed around 2 K. 
The next section is devoted to the analysis of the suscep- 
tibility in the paramagnetic regime above the anomaly, 
the nature of the phase below it will be the subject of 
another publication [23j . 

The inverse x~^ of the initial magnetic susceptibility, 
shown in Fig. \5\ has a peculiar shape, where three dis- 
tinct regions, with different (almost) linear slopes at low 
(I), intermediate (II), and high (III) temperatures, can be 
distinguished. In a first analysis following ref. [2l|, re- 
gions (I) and (III) were assumed to correspond to two dis- 
tinct paramagnetic regimes of interacting magnetic enti- 
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FIG. 5: La3Cu2V09 inverse linear magnetic susceptibility 
measured in 0.1 T and 6 T below and above 350 K respectively 
(open circles). Three regions of distinct magnetic behaviours 
are materialized in grey. A low temperature zoom is shown in 
the inset. The straight lines in regions I and III are calculated 
in a Curie- Weiss model (see text). 



Pulsed fields 




H(T) 

FIG. 6: Magnetization isotherm measured around T=1.4 K 
in increasing and decreasing pulsed field up to 55 T, and its 
field derivative. The dotted line indicates the middle of the 
magnetization plateau, which is related to the Zeeman cross- 
ing of the two first levels of the system with AS=-|-1 (see inset 
of Fig. 0. 



ties accounted for, in the mean field approximation, by a 
Curie- Weiss law x = C/{T—6) with the Curie- Weiss tem- 
perature 6. The Curie constant, C = ng^fiBS{S+l)/3kB 
with the Lande factor g=2, is evaluated for 1/3 of 8-spin 
cluster -I- 2/3 of 9-spin clusters, i.e. 8.667 Cu^+; n is 
therefore the number of S=l/2 spins per cluster. The 
resulting Curie- Weiss analysis are shown in Fig. [5] and 
the numerical values are reported in Table |TT| n, 6, the 
effective moment ^^ff — g\/S{S + 1), and the temper- 
ature range of the best fits corresponding to the linear 
portion of regions (I) and (III). The slope of the high 
temperature region (III) is well described above 500 K 
with a Curie- Weiss fit using the expected value n of 8.667 
5=1/2 spins per cluster in the paramagnetic regime (Fig. 
[5]) . The resulting large negative Curie- Weiss temperature 
for this regime, 9 = —294 K, denotes strong intra-cluster 
antiferromagnetic couplings between these spins. 

In the low temperature region (I), it is more difficult 
to isolate a linear regime in 1/x, which rather presents 
a continuous curvature (Fig. |4]). A forced Curie- Weiss 
fit in a reduced temperature range yields approximately 
one S—1/2 spin per cluster. The reduction of the Curie 
constant, by a factor close to 9, suggests that the mag- 
netic entities in the low temperature range are collective 
5=1/2 pseudo-spins, resulting from the entanglement of 
the wave functions of the paramagnetic spins within each 
cluster at high temperature. These pseudo-spins seem to 
be weakly antiferromagnetically coupled to each other 
since 9 = —4.2 K in this region, which is about two 
orders of magnitude weaker than the 9 value of region 
(III) corresponding to the intra-cluster coupling. This 
value, close to the one at which the susceptibility presents 
an anomaly, suggests that this last feature is related to 
the inter-cluster coupling which becomes effective at this 
temperature (cf . section IIIB II) . We will come back to 



FIG. 7: High field magnetization isotherms in meaningful 
units, measured by extraction at 3 K , and compared to the 
BS magnetization measurements. The solid and dotted hor- 
izontal lines point out the calculated saturation value of the 
ground state and of the first excited state magnetizations (see 
section HIT)) . The Zeeman diagram corresponding to the field- 
induced crossing of these two levels is shown in the inset. 



the analysis of the low temperature regime in section IIIII 
using a more detailed description of the clusters. 

Concerning the intermediate temperature region (II), 
although 1/x could also be believed to vary linearly from 
40 to 200 K [HI, a Curie- Weiss analysis does not have 
any physical meaning since no collective magnetic entities 
are formed, as will be shown in the theoretical analysis 
of section IIIII 



6 



Region T range (K) 9 (K) n^ff (i^b) n 
I 2-11 -4.2 ± 1.1 1.732 1.02 ±0.19 

III 500-850 -294 ± 73 1.732 8.667 



TABLE II: Curie- Weiss parameters of the best fits accounting 
for the inverse susceptibility in the two quasi-linear regimes 
at low and high temperature. In the high temperature region, 
the number of spins 5'= 1/2 per cluster was hold constant and 
equal to its expected value of 8.667. 



2. High field magnetization 

High field magnetization measurements were per- 
formed using pulsed fields at the Laboratoire National 
des Champs Magnetiques Pulses (Toulouse). The pulsed 
field, which can reach 55 T, is produced by the discharge 
of a capacitor bank into a resistive copper coil. The 
signal in the presence of the sample is proportional to 
the time derivative of its magnetization. The M vs H 
curve obtained at T « 1.4 K is displayed in Fig. [S] 
It shows a saturation plateau centered on the inflexion 
point Hq/2 ss 29 T, indicated by the minimum in the 
magnetization field derivative. Note that the influence of 
the weak inter-cluster interactions should be irrelevant 
at such high magnetic fields. An upturn of the magneti- 
zation occurs at higher magnetic fields, suggesting that 
the plateau should end at about Hq = 58 T. This is 
attributed to a Zeeman crossing of the first excited mag- 
netic energy levels with the ground level of the cluster 
system (cf. inset of Fig. [7]). The quantitative analysis of 
this crossing will be further discussed in section IIIII 

Since the pulsed field technique does not allow abso- 
lute measurements with a good accuracy, complementary 
measurements at 3 K were performed at the Grenoble 
High Magnetic Field Laboratory (GHMFL) on a magne- 
tometer using the axial extraction under a magnetic field 
up to 23 T produced by a 10 MW resistive magnet. An 
absolute calibration was performed using a Ni sample. 
The resulting curves presented in Fig. [71 in agreement 
with the BS lower fields measurements, show a tendency 
to saturate at a value close to 0.67 ^s/cluster, the ex- 
pected value for antiferromagnetic 8- and 9-spin clusters 
in the (l/3)/(2/3) ratio (cf. section . 



3. Specific heat 

Specific heat Cp in zero magnetic field was measured at 
the CEA-Grenoble on a commercial calorimeter (Quan- 
tum Design PPMS) using the pulsed relaxation method 
and equipped with an '^He refrigerator allowing to reach 
400 mK. Its low temperature behavior is shown in Fig. 
[HI An anomaly is clearly visible at 1.82 K, coinciding 
with the one observed in the initial magnetic suscepti- 
bility X- In the absence of a non-magnetic isostructural 
compound, a determination of the magnetic part of the 
specific heat in the full temperature range is difficult. 




T(K) 

FIG. 8: (Color online) Specific heat measurements at low tem- 
perature (circles) analyzed with two models: (i) The phonon 
contribution (dashed line) is obtained from a fit of Cp/T as 
a function of using the low temperature approximation 
Cp = /JiT^ + l32T^ in the range [12-28 K]. The resuhing mag- 
netic signal is represented with square symbols, (ii) The sec- 
ond model is based on the multi-J spin cluster description 
reported in table IIIII of section IIIII The lattice contribution 
of the form Cp = (3iT^ + f32T^ (dotted lines) is adjusted such 
that the remaining signal (triangles) coincides with the cal- 
culated magnetic specific heat (solid line) above 5 K. 



However, the lattice contribution is often evaluated by 
a fit of Cp/T versus in a temperature range where 
it is assumed that the magnetic signal is negligible and 
that the following low temperature approximation for the 
phonon contribution is valid: Cp = (3iT^ + /32T^ with 

f3i = ^^^NkB/Q% where Qd is the Debye temperature 
and N is the number of vibrating ions. In the present 
case, the best fit in the 12-28 K temperature range yields 
the values: f3i = (7.3 ± 0.09)10^'*, as from which we 
get Bd ~ 271 K, and (32 = (-2.3 ± 0.1)10-^. The as- 
sumed magnetic part of the specific heat, remaining after 
subtraction of this contribution, presents a broad signal 
centered around 6 K, as shown in Fig. [51 This signal is 
interpreted as the signature of a Schottky anomaly due 
to the presence of discrete low energy levels associated to 
the spin clusters. This analysis will be deepened using 
the finer description of the pseudo-spins obtained from 
the calculation discussed in section IIIII 



4-. Inelastic neutron scattering 

In order to probe the energy spectrum of the assembly 
of collective spins in La3Gu2V09 , inelastic neutron scat- 
tering measurements were performed on the IN4 time-of- 
fiight spectrometer at the Institut-Laue-Langevin with 
an incident neutron beam of energy 16.9 meV (i.e. of 
wavelength A=2.2 A), and an energy resolution of 0.77 
meV. The spectra were corrected from the background 
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contribution and the detector efficiency was calibrated 
using a vanadium sample. 

The energy response recorded at 2, 150 and 300 K 
integrated over the small momentum transfer Q, more 
precisely in the range [0.1, 2] A"-'^, where the magnetic 
scattering is expected to be the largest are displayed 
in Fig. [9l As apparent from the increase of the sig- 
nal with increasing temperature and with increasing Q, 
there is an important contribution from the scattering by 
the phonons in the investigated energy range (up to 15 
mcV on the neutron loss energy side). A way of remov- 
ing it is to substract the spectrum at high temperature 
T//=300 K renormalized by the appropriate thermal fac- 
tor (l-exp(--^|f^))/(l-exp(--^|f^)) from the spectra 
at lower temperatures T/^. This treatment might under- 
estimate the magnetic contribution since magnetic signal 
arising from transitions between energy levels can still 
be present at high temperature, especially if the energy 
spectrum of the system is dense and extended in energy. 
This analysis allowed nevertheless to visualize in the cor- 
rected spectra a broad bump around 8.7 meV at 2 K , 
which largely decreases at 150 K. The magnetic origin 
of this feature was confirmed by its slow decrease with 
increasing Q from 1 to 4 A~^, once corrected from the 
phonons. An additional signal at 2 K seems to be present 
in the foot of the elastic peak. The nature of this signal 
is hard to establish from this sole measurement since it 
could as well be quasielastic or inelastic with an energy 
gap of the order of the energy resolution. 
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FIG. 9: (Color online) Inelastic neutron scattering measure- 
ments. Upper frame: Energy spectra at 2, 150 and 300 K, 
integrated on the Q interval [0.1, 2] A~^. Lower frame: 2 and 
150 K spectra corrected from the phonon contribution using 
the 300 K spectrum (see text). 



III. ANALYSIS AND DISCUSSION 

The crystallographic study established that each 
La3Cu2V09 plane can be described as a triangular lat- 
tice of randomly distributed 1/3 of 8-spin and 2/3 of 
9-spin clusters with antiferromagnctically interacting lo- 
calized S'=l/2 spins on each Cu. The magnetic proper- 
ties of La3Cu2V09 in the range of temperature T > 2 
K should be dominated by a physics of isolated clusters, 
since the inter-cluster interactions are two orders of mag- 
nitude smaller than the intra-cluster ones. In the follow- 
ing, we shall probe this regime through calculations by 
exact diagonalization of 8-spin and 9-spin clusters model. 
The intra-cluster interactions parameterizing these calcu- 
lations will be adjusted from comparison with the exper- 
imental results. These impose the following constraints: 
first, the calculation should reproduce the measured ther- 
mal variation of the initial magnetic susceptibility and 
the magnetization processes in high magnetic field, in- 
cluding the saturation plateau and the Zeeman crossing 
at the largest available magnetic field. Second, the cal- 
culated energy spectrum of the system should account 
for the energy excitations revealed by inelastic neutron 
scattering and the magnetic specific heat. The nature of 
the resulting pseudo-spin state at low temperature will 
then be discussed. 



3500-1 




T(K) 

FIG. 10: Comparison of the inverse of the measured suscep- 
tibility (circles) with calculated ones (line) from a single J 
model (grey) and from the multi-J model of table Hill taking 
into account 8-spin (dotted) and 9-spin clusters (dashed). 



A. Cluster calculation versus experiments 

The calculation of the La3Cu2V09 magnetic properties 
are obtained through exact diagonalisation of the Hamil- 
tonian of the model system where each 8-spin and 9-spin 
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FIG. 11; Calculated total M(H) in the multi- J model of table 
Hm with the 8 -spin (dotted) and 9-spin (dashed) contribu- 
tions, at 5 K. The grey line is the field derivative of the total 
M(H). 



cluster is described by the Heisenberg Hamiltonien 



(1) 



where the Jij < 0, labeled in Fig. [31 stand for all near- 
est neighbor antiferromagnetic interactions between the 
spins i and j of the cluster (Fig. [H]). In a single- J model 
assuming that all Jij = J = —445 K, the ground states of 
the 9- and 8-spin clusters are found to be associated with 
a total spin S—1/2 and S=0, respectively. This yields in 
average 0.667 5=1/2 spin per cluster at low temperature, 
which is a bit smaller than the value extracted from the 
Curie- Weiss fit of the paramagnetic regime at low tem- 
perature (dotted hne in the inset of figure [S]). The mag- 
netic susceptibility and other thermodynamics properties 
of the system were computed and compared to the ex- 
perimental data. A single-J model does not account for 
the initial magnetic susceptibility as shown in Fig. \W\ 
where one can see strong deviations from the experimen- 
tal data below T « 300 K down to low temperatures. 
This demonstrates the need to consider different J val- 
ues, that will be expressed, in the following, in units of 
J = -445 K. 

To push further the analysis, some additional assump- 
tions are required in order to reduce the large number 
of non- equivalent first-neighbor interactions : 4 and 8 
distinct Jij values for the 9- and 8-spin clusters respec- 
tively, as shown in Fig. [31 In the 8-spin clusters, the 
main consequence of the distortion induced by the Cu/V 
substitution is an elongation of the J7 and Js exchange 
paths (see section III A 21 and Fig. [3]). These exchange 
parameters are thus expected to be much weaker than 
all the others in the 8- and 9-spin clusters. We therefore 
assume a 3-J model for the 8-spin clusters, i.e a model de- 
fined within the ( Ji — . . . — Jg, J7, Jg) parameter space. 
Concerning the 9-spin clusters, there are no straightfor- 



TABLE III: Example of multi- J model that accounts well 
for all the experimental data. The 8-spin and 9-spin clusters 
nearest-neighbor interactions J, labeled as in figure [S] are 
given in units of 445 K. 



ward assumptions that could be made and the 4-J model 
described in Fig. [3lis adopted. 

The measured magnetization processes under magnetic 
fields provide with meaningful constraints over the J7 and 
Js parameters. The saturation of the magnetization, cor- 
responding to the alignment of the S—1/2 pseudo-spin of 
the 9-spin clusters anti parallel to the applied magnetic 
field, is reproduced already within the single J model. As- 
suming weaker J7 and Jg, the first excited states which 
will be field-driven to energy lower than that of the zero 
field ground state are those belonging to the 8-spin clus- 
ter excited manifold of states with total spin S—l and 
describing antiparallel alignment to the applied magnetic 
field. The corresponding Zeeman crossing experimentally 
inferred to occur at about 58 T then imposes that J7 and 
Jg cannot take values smaller than 0.28 and 0.16, respec- 
tively. In order to get the inflexion point Ho/2 around 
29 T, one of these parameters must at least take the 
above limit values, as shown in Fig. [TTlfor J7=0.28 and 
Jg=0.24. Note that M(II) has been calculated for T=5 
K in order to round the sharp edges of the calculated 
magnetization saturation plateau and thus better repro- 
duce the measurements. This shape difference together 
with the fact that the experimental saturation is slightly 
higher than the calculated one could be due to a distri- 
bution of J values not taken into account in our simplest 
model (some kind of disorder, influence of inter-cluster 
coupling, presence of a tiny amount of clusters with dif- 
ferent Cu/V substitutions etc.) 

The excitation spectrum associated with this 8-spin 
cluster model is shown in figure [121 A first bunch of 
levels is observed around 8.5 meV above the ground state, 
and the next bunch lies at energies greater than 15 meV. 
At 2 K, this first bunch of levels will therefore yield an 
excitation at an energy value in agreement with the one 
observed in inelastic neutron scattering experiment (cf. 
Fig. ^. 

Small variation of the Js on the other bonds of the 
8-spin and 9-spin clusters are expected from the slight 
differences in their Cu-0 distances and Cu-O-Cu angles 
and by the different coordination environments of the 
Cu(2) site on one hand and of the Cu(3) and Cu(4) sites 
on the other hand (section [II A 1[) . An experimental con- 
straint that gives some information about this distribu- 
tion of J values is the broad bump around 6 K observed 
in the magnetic specific heat. No magnetic contribution 
to the specific heat is observed below 10 K in the single-J 
model. Conversely, small differences, no larger than 5 %, 
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FIG. 12: 8-spin and 9-spin clusters degeneracies as a function 
of the energy difference with respect to the fundamental level 
(Eo) in the low energy region, for the multi-J model of table 
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between the three bonds of the external triangles in the 
9-spin clusters allow the rise of a bump in the calculated 
specific heat around 6 K. Note that the low energy spec- 
trum as well as the shape of the magnetic susceptibility 
are much less sensitive to the value of J4 (the inner tri- 
angle exchange) . In addition, similar variations of the Ji 
to Jq bonds in the 8-spin clusters does not yield any fea- 
ture in the specific heat in this low temperature range. 
On this basis, a multi-J model with Ji=1.0, J2=1.01, 
J3=1.045, and J4=1.02 in the 9-spin clusters (cf. Ta- 
ble |TTT| inferred from the specific heat, was proposed and 
found to also well reproduce the inverse magnetic suscep- 
tibility, as shown in figure \W\ 

The resulting specific heat is shown in figure [H] In ad- 
dition to the requested bump around 6 K, the calculation 
also reveals the presence of significant magnetic signal at 
high temperatures. This implies that the evaluation of 
the phonon contribution that was made section III B 3[ 
assuming that there is a region where it follows the low 
temperature Debye approximation and where the mag- 
netic contribution is negligible, is incorrect. However, 
a phonon contribution of the form Cp — (3iT^ + P2T^ 
with Pi = 0.00047, corresponding to 8_d « 314 K, and 
/32=1.1184 10~^ can be adjusted such that the remain- 
ing signal coincides with the calculated magnetic specific 
heat above 5 K (cf. figure [5]). 

The low lying energy levels obtained through such 
small variations of the exchange parameters on the ex- 
ternal triangle in the 9-spin clusters are shown via the 
histograms of figure [T^ They would also be in agreement 
with the presence, in the neutron scattering experiment, 
of a magnetic inelastic signal at low temperature in the 
foot of the elastic peak (cf. figure ^ . 

Although not necessary for the matching of the calcula- 
tions with the experimental data, small variations of the 
Ji to Jg exchange parameters in the 8-spin clusters are 



also expected to be present in the real system. The exact 
distribution of all the 8-spin and 9-spin clusters exchange 
parameters can however not be accessed unambiguously 
from the present data. Because of the large number of 
adjustable parameters, the multi-J model presented ta- 
ble lllll is thus not the only relevant, but the main features 
of this model, arising from the experimental constraints, 
are quite robust : smaller values of J7 and Jg in the 
8-spin clusters with respect to all other first-neighbors 
parameters, and distribution of exchange parameter val- 
ues of the order of less than 5 % in the 9-spin clusters. 
Note that this model gives a satisfactory agreement with 
all the experimental data. The inverse magnetic suscep- 
tibility in particular is much better reproduced than in 
the single- J model, except for tiny differences at low tem- 
peratures. These last discrepancies could be due to the 
onset of inter-cluster correlations at low T that were not 
included in the model. 

The evolution of these pseudo-spins when lowering 
temperature can be probed through the square of the 
total spin of the cluster (S*^). Commuting with the clus- 
ter Hamiltonian, (S'^) is a good quantum number that 
allows to track the spin value of the entangled state in 
the whole investigated temperature range. In the multi-J 
model, the pseudo-spins are stabilized below w 18 K for 
the 9-spin clusters and below « 8 K for the 8-spin clusters 
as shown by the constant (5^) values, 3/4 and respec- 
tively, observed below these temperatures (cf. inset of 
Fig. I13p . Concerning behavior at higher temperatures, a 
1/T expansion of the calculation in the region III, con- 
firms that the system can be described by a Curie- Weiss 
analysis in the temperature range T ^ 550 K as done 
in section Tub 11 In the intermediate temperature region 
II, no plateau regime interpretable in terms of entangled 
spins states nor a 1/T behavior is observed. The vari- 
ation of the magnetic susceptibility in this temperature 
regime (cf. figure [5]) is rather related to the onset of the 
strong magnetic intra-cluster correlations and cannot be 
ascribed to any building of pseudo spins. 



IV. CONCLUSION 

The La3Cu2V09 oxide compound is shown to be con- 
stituted of 8-spin and 9-spin clusters of 5'= 1/2 spins, laid 
out on 4 vertex-sharing triangles, building block of the 
kagome lattice. From susceptibility, high field magnetiza- 
tion, specific heat and inelastic neutron scattering mea- 
surements analyzed with exact calculation of spin cluster 
models, the low temperature stabilization of collective 
pseudo-spins S'=0 and S=l/2 in each 8-spin and 9-spin 
cluster, is revealed. From Susceptibility and specific heat 
measurements, the coupling of these collective entities 
is evidenced below 2 K. This study opens the way to 
original measurements of the individual and cooperative 
behaviour of collective quantum entities, for instance a 
direct imaging of their wave function using polarized neu- 
tron scattering. 
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